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ABSTRACT 

We constrain the distance of the Gamma-Ray Burst (GRB) prompt emission site from the 
explosion centre, R, by determining the location of the electron's self absorption frequency 
in the GRB prompt optical-to-X/7-ray spectral energy distribution, assuming that the optical 
and the 7-ray emissions are among the same synchrotron radiation continuum of a group of 
hot electrons. All possible spectral regimes are considered in our analysis. The method has 
only two assumed parameters, namely, the bulk Lorentz factor of the emitting source V, and 
the magnetic field strength B in the emission region (with a weak dependence). We identify 
a small sample of 4 bursts that satisfy the following three criteria: (1) they all have simulta- 
neous optical and 7-ray detections in multiple observational time intervals; (2) they all show 
temporal correlations between the optical and 7-ray light curves; and (3) the optical emission 
is consistent with belonging to the same spectral component as the 7-ray emission. For all the 
time intervals of these 4 bursts, it is inferred that R ^ 10 14 (r/300) 3 / 4 (B/10 5 G) 1 / 4 cm. 
For a small fraction of the sample, the constraint can be pinned down to R rj 10 14 — 10 15 
cm for r ~ 300. For a second sample of bursts with prompt optical non-detections, only up- 
per limits on R can be obtained. We find no inconsistency between the i?-constraints for this 
non-detection sample and those for the detection sample. 

Key words: gamma-rays: bursts - gamma-rays: theory - radiation mechanisms: non-thermal 
- radiative transfer 



1 INTRODUCTION 



Although the Gamma-Ray Burst (GRB) was discovered about 50 
years ago first through its prompt 7-ray emission, large uncertain- 
ties still remain in understanding the prompt emission site, namely, 
the distance of the emission site from the explosion centre R, with 
controversial evidence. There are three possible sites discussed in 
the literature. One is the standard internal- shock site which depends 
on the fluctuation time scale 5t seen in GRB light curves (e.g., Rees 
& Meszaros 1994, see Piran 2005, Meszaros 2006 for reviews). It 
can have a large range of R ~ T 2 c8t ~ 10 13 - 10 15 cm because 
St and r vary largely from burst to burst. The second is the pho- 
tospheric radius at 10 11 — 10 12 cm at which the prompt emission 
arises as a combination of the photosphere thermal emission and 
a Comptonized component above it, the latter being induced by 
some energy dissipation process below and above the photosphere 
(e.g. Rees & Meszaros 2005; Ryde et al. 2006; Thompson et al. 
2007). The third one is a large radius (> 10 14 cm) as is supported 
by the Swift XRT data (Lazzati & Begelman 2005; Lyutikov 2006; 
Kumar et al. 2007) and Fermi data of GRB 080916C (Abdo et al. 



* E-mail: rfshen@astro.as.utexas.edu (R-FS); zhang@physics.unlv.edu 
(BZ) 



2009; Zhang & Pe'er 2009), possibly due to magnetic dissipation 
(e.g., Lyutikov & Blandford 2003). 

The rapidly responding ability of a few GRB-dedicated 
ground or space based optical telescopes, e.g., ROTSE (Akerlof et 
al. 2003), RAPTOR (Vestrand et al. 2002), TORTORA (Racusin et 
al. 2008) and the UVOT (Roming et al. 2005) on aboard the Swift 
satellite, has enabled the time-resolved detection of bright prompt 
optical emission before the 7-rays die off, for about a dozen of 
GRBs. Five of these GRBs, i.e., 041219A (Vestrand et al. 2005), 
050820A (Vestrand et al. 2006), 05 1 1 1 1 (Yost et al. 2007a), 061 121 
(Page et al. 2007) and 0803 19B (Racusin et al. 2008), show a tem- 
poral correlation between the strongly variable optical flux and the 
7-ray pulses, which suggests that the optical emission most likely 
shares the same dynamical process that is responsible for the highly 
variable 7-ray emission. While the other four bursts have optical 
flux densities below or marginally consistent with the extrapola- 
tions from the low-energy power law of the 7-ray spectra, the op- 
tical flux density in GRB 0803 19B exceeds the 7-ray extrapolation 
by 4 orders of magnitude (Racusin et al. 2008; Kumar & Panaitescu 
2008), suggesting that for this burst alone the optical emission has 
a spectral origin different from that of the 7-rays. 

In this paper, for the four GRBs - 041219A, 050820A, 051 1 1 1 
and 061121 - we assume that the prompt optical and the 7-ray 
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emissions are components belonging to the same synchrotron radi- 
ation continuum of a group of hot electrons. Based on this assump- 
tion, the self-absorption frequency of the synchrotron electrons, u a , 
which causes a break in the long-wavelength part of the continuum, 
can be determined or constrained by studying the optical-to-7-ray 
spectral energy distribution (SEDjj. Since u a is dependent on the 
properties of the prompt emission source, such as the distance of 
the emission site from the explosion centre R, the bulk Lorentz fac- 
tor (LF) r and the magnetic field B of the source, from u a we can 
determine or make constraints on R for these bursts, using infor- 
mation on r and B obtained in other ways. This is the main goal of 
this paper. Since the prompt optical and 7-ray components in GRB 
0803 19B are most likely of different spectral origins because of its 
peculiar SED shape, our approach is not applicable to this burst. 

On the other hand, for some other long GRBs the rapid re- 
sponse of the dedicated ROTSE telescope has returned only upper 
limits of the optical flux density during the prompt phase (Yost et 
al. 2007b). Another goal of this paper is to get constraints on R for 
these optically "dark" bursts and to study whether the prompt op- 
tical non-detection is caused by a heavier self-absorption due to a 
closer emission site to the explosion centre. 

In this paper, we first derive analytically u a in terms of R, T, 
B and the emission properties in Sec. 2. The arguments that sup- 
port our assumption of one synchrotron continuum component for 
both optical and 7/X-ray are given in Sec. 3. We derive in Sec. 4 
the constraints on R through u a explicitly, by determining the loca- 
tion of v a in the optical-to-7-ray SED and considering all possible 
spectral regimes. We apply this method to a prompt optical detec- 
tion GRB sample and a prompt optical non-detection sample which 
are described in Sec. 5. The results are presented in Sec. 6. Finally 
the conclusion and discussions are given in Sec. 7. 



2 DETERMINING THE SELF-ABSORPTION 
FREQUENCY 

The GRB high energy emission spectrum is characterized by a 
smoothly joint broken-power-law form (Band et al. 1993). Thus the 
relativistic electrons responsible for the GRB prompt emission due 
to either synchrotron or synchrotron self-inverse Compton (SSC) 
radiation are in a piece-wise two-power-law energy distribution: 



7V(7) oc 



if 7m < 7 < 7p • 
if 7 > 7p , 



(1) 



where N(y) is such defined that N^dy is the number density of 
electrons with energy in the interval of 7 to (7 + d'y), and y m is the 
minimum energy of these relativistic electrons (for convenience we 
omit the factor m e c 2 in the electron energy 7m e c 2 throughout the 
text when electron energy is mentioned). 

Note that although this distribution set-up is phenomenologi- 
cally based on the two-power-law shape of the high energy radia- 
tion spectrum observed in GRBs, it has specific physical meanings. 
Within the shock acceleration scenario, newly accelerated electrons 
with a minimum energy 7, and a power-law energy distribution 
are continuously injected. These electrons lose energy through ra- 
diative cooling, and the instantaneous electron spectrum steepens 
above a critical energy 7 C . All the electrons with energy larger than 



1 The significance of self-absorption frequency has been highlighted by 
Doi, Takami & Yamazaki (2007) who used the varying location of u a to in- 
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7 C radiate away their energy within a time shorter than the dynam- 
ical time. When 7,: < y c , our notation corresponds to y m = 7, 
and 7 P = y c . When y c < 74, the cooling causes a flatter power 
law between y c and 7^, even though the newly accelerated elec- 
trons are injected in the energy range above 7*. For this case one 
has 7 m = 7 C and 7 P = 7i. In summary within the shock acceler- 
ation scenario one has y m = min(7i, y c ) and y p = max(7i, y c ). 
More generally, one can also have a scenario that invokes contin- 
uous heating and cooling of electrons (e.g. that envisaged in the 
reconnection models), and y p then reflects the intrinsic break in the 
steady state electron spectrum. In any case, our treatment is generic, 
which does not depend on the concrete particle acceleration mech- 
anism and the origin of 7 P . 

The broken power-law electron energy spectrum naturally 
gives rise to a piece- wise power law photon spectrum as commonly 
observed: 



f v oc 



if u m < v < Up 
if u > u„ , 



(2) 



where /„ is the observed flux density (in units of mjy), u m and 
u p are the observed characteristic emission frequencies of the elec- 
trons with energy j m and j p , respectively, and u p is usually the 
peak frequency of the uj v spectrum. The low-energy power law 
u 131 does not extend to low frequency indefinitely. Without syn- 
chrotron self absorption, the spectral index below u m changes to 
1/3, regardless of whether 7^ < 7 C or 7 C < ft. Below a certain 
frequency u a <C u p , the synchrotron self absorption starts to play a 
significant role - at frequencies lower than u a the emitted photons 
are thermalized with the electrons. The self-absorption frequency 
u a is such defined that at this frequency the self-absorption optical 
depth T sa {u a ) = 1. 

Let us determine u a for an emitting GRB ejecta moving with a 
Lorentz factor (LF) T at a distance R from the center of the explo- 
sion. In the ejecta comoving frame (hereafter the quantities mea- 
sured in this frame are marked with the prime sign), u' a can be de- 
termined by equating the un-absorbed source surface flux density, 
F' v , , at u'a to a blackbody surface flux density with temperature T' 
in the Rayleigh- Jeans regime (see Appendix for the derivation; also 
see Sari & Piran 1999, Li & Song 2004, McMahon, Kumar & Piran 
2006): 



2kT 



C(Pi)F', 



for v' m < u'a < u' p , 



(3) 



where T 1 = max(7 a , 7 m )m e c 2 /A:, y a is the energy of electron 
whose characteristic emission frequency is u a , and k is the Boltz- 
mann constant. The numerical factors Ci and C2 are functions of 
f3i only whose values range from 1.2 to 4.5 and from 1.2 to 7.0, re- 
spectively, for the range of observed /3i values, but they have been 
neglected in previous works while we include them here. 

Transforming the frequency to that measured in the observer's 
frame gives u a = u' a T /{I + z), where z is the redshift of the GRB 
host galaxy. Measuring in the host comoving frame, the source 
has an isotropic luminosity of 4irR 2 T 2 F', uL. This luminosity is 

also given by A-kD\J Vv u p , where Dl is the luminosity distance of 
the GRB and f v is the observed peak flux density. Thus we have 
K> = fv P ( D L/R) 2 /[F(l + z)]. After applying these relations, 
Eq. (3) becomes 
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for f a < I> m ! 



+ z) J 



2 



1 



m e 7 i/^ (-^-) 

for V m < Va 



R \2 (1 + z) 3 



(4) 



After substituting the electron's energy 7 using the following 
relation between 7 and the photon frequency v, 



1 = 



) 2 (n^) 2 for synchrotron, 

1 i/i , for SSC, 



' 27rm e c ^ 2 

„ eB 

/ 27rm e c \ 4 / l + z \ 4 



(5) 



k V eB } \ r 
the self-absorption frequency is calculated as: for synchrotron 

fa = 

l 2 J X 1U /"P^p,19 "m,17 



38. 5-190! 



Hz 

if f a < v m ; 
-Pi 



(6) 



28 s ) 2.5-01 / r 3 QD \r5^7g5-2ft ^2.5-01 



1 + 2 . 



Hz 



if f m < V a < Vp 



and, for SSC 



'2l)« x 10 15 -t^/I^o /3l K sA 



• D 



L .28 



1 + z / V 1 + 

35-1901 



1, 

5 A 

p"p, 19 "m,17 

f(r^)l B A^-f Hz 

if U a < V m \ 



(C2.y71S=Wi x 1Q 2.25-0! f£- , ,, 



-01 



2.25-0! 2.25-01 
"p,19 

°L,28 ~\ 2.25-01 / T l 300 \ 2"7#^37 R 9-401 p2. 25-01 
1 + z J I 1 + z J -°5 rt 14 



(7) 



K B^ P1 R^ a - pi Hz 

if V m < V a < Vp ■ 



In the results above, the flux density, e.g., /„ , is in units of mjy, 
T = 300xT3oo and the convention Q = Q n X 10™, e.g., v = v 19 x 
10 19 Hz and B — B$ x 10 5 Gauss, is used for other quantities; the 
same notations will be used in the rest of the paper. In the following, 
our discussion will be based on the synchrotron radiation only. But 
for the use of reference the expression of v a for SSC is also given 
here. 



3 THE ONE-SPECTRAL-COMPONENT ASSUMPTION 

In this paper, we make an assumption that for the two samples (see 
Sec. 5 for a description of the sample selection criteria) studied, 
the optical and 7/X-ray photons belong to a same synchrotron con- 
tinuum spectrum generated by a same group of hot electrons. This 
assumption is based on the following three considerations. First, 
GRB prompt 7/X-ray emission is often interpreted as synchrotron 
radiation of a group of non-thermal electrons. If that is the case, the 
synchrotron spectrum must have a continuum extending to the low- 
frequency regime up to the optical band, presumably with a gentle 
slope of f v oc v 1//3 unless it has a self-absorption break. Secondly, 
the optical flux density is expected to always lie below or near the 
extrapolation from the 7-ray spectrum. This is generally consistent 
with observations in our sample (see Sec. 5 below and Fig. 2). Fi- 
nally, a temporal correlation between the optical flux variation and 
the 7-ray LC is observed for the 4 GRBs in the our sample, sug- 
gesting the two components likely have the same dynamical origin. 
This is the major supporting evidence for our assumption. 



We notice that there are other scenarios on prompt 7/optical 
emission that have been discussed in the literature. These include 
the synchrotron + SSC model (Kumar & Panaitescu 2008; Racusin 
et al. 2008), the models invoke different emission radii for opti- 
cal and 7-ray emissions (Li & Waxman 2008; Fan, Zhang & Wei 
2009), and the model invokes two shock regions at a same emis- 
sion radius (Yu, Wang & Dai 2009). These models are relevant to 
GRB 0803 19B, which clearly requires a separate spectral compo- 
nent to interpret the prompt optical emission. For most other bursts 
studied in this paper, although these models are not ruled out, they 
are not demanded by the data. Our simple one-component model is 
adequate to interpret these bursts, and we will hereafter adopt this 
one-component assumption. 



4 DERIVING R CONSTRAINTS FROM v A AND SED 

Let us consider the synchrotron emission as the prompt emission 
mechanism. Depending on the locations of v m and v a , the ratio of 
the optical flux density to the flux density at v v is given for four 
different spectral cases by 

fvopt _ 



Pi 



if v m < v pi, < v a 

if V m <V a < V opt 



-) 3 , if V a < V opt < V, 



< V a < V „ 



case (I), 
case (II), 
case (III), 
case (IV). 



Fig. Q] illustrates these four spectral cases. There is also 
a variation of case (III) (let us call it case III. 5): v opt < 
v m < v a , for which the flux density ratio is fu opt /fv p = 
(va/vp) 131 (v m /v a ) 5 / 2 (v pt/v m ) 2 . We will come back to this case 
later and show that it is almost exactly the same as case (III). In 
some GRBs simultaneous observations of prompt X-ray and 7-ray 
emissions are made. The X-ray spectrum corrected for the photo- 
absorption in the soft end always nicely matches the power law 
extrapolated from 7-ray spectrum below v p , without a need to in- 
voke a break (e.g. Cenko et al. 2006; Romano et al. 2006; Page et 
al. 2007). This requires that both v a and v m be smaller than vx = 
0.3 keV, the lower end of the Swift XRT band pass. We take this as 
a constraint in our analyses. 

We aim to constrain R based on the spectral information, such 
as f Vp , f v t and f3i, using v a as a proxy. We know from Sec. 2 
that v a is expressed in terms of f Up , v m and R. So for each spec- 
tral case, we substitute the appropriate v a expression into the flux 
density ratio equation or into the constraint on v a implied by the 
definition of that spectral case (Eq. 8), and then get the R con- 
straints. 

For case (I), substituting v a and letting v op t — 5 x 10 14 Hz 
(for R band), we have 

« U =7.5x(f)^i,^4(^)». (9) 

There is also a justification criterion due to the case definition v a > 
v pt, where v a is directly determined from the flux density ratio in 
Eq. (8) by 



10 



(19- 



■2 r, 1 



p,19 



fv p 

fvopt 



Hz. 



(10) 
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Log v Log v 



Figure 1. The four cases of the broad-band synchrotron spectrum of the GRB prompt emission discussed in the text. The broken power law shape around v p 
is phenomenologically derived from the observed 7-ray spectrum in GRBs; see its possible theoretical origins in Sec. 2. u a is the self absorption frequency. 
The dotted line in the two top panels is the low energy part of the spectrum when the self absorption is not considered. The two grey vertical bars marks the 
positions of v pt and vx, respectively, ux = 0.3 keV is the lower end of the Swift XRT band pass. 



For case (II), the flux density ratio does not depend on v a and 
hence on R. One justification criterion for this case is that the spec- 
tral slope from the optical to the X- or 7-rays has to be consistent 
with f3\, i.e., /3 opt _x/ 7 = /9i. Another criterion due to the case 
definition is v a < v ovt . Substituting with the expression of v a , the 
latter gives 

fl 14 >7.5x (Q.y x(2.2) dl xl(T 2 - 5/3 i 

Dl,28 



ih, 

2 

'p,19 ■ 



I J \l + zj 



(11) 



For case (III), when substituting the appropriate v a expression 
in Eq. (6) into the flux density ratio relation, it gives the expression 
of R in 



Ri 



2-0x(-j u^f^^-B^—) . (12) 

We find for case (III.5) that the R expression — obtained by sub- 
stituting the appropriate v a expression into Eq.(8) — is exactly the 
same as Eq. (12) except that C\ is replaced with C2. We find the 
ratio C2/C1 lies in the range of (1, 1.6) for /3i — ( — 1.4, 0) (see in 
Appendix). 

According to the definition of case (III), v op t < v m < vx- 
Plugging this constraint of v m into Eq. (12) and its counterpart 
equation for case (III. 5), it gives 

2 - lx (-) f ^—) <R ^ i {T^) < 



'■ bX \2) U ^ (1 - 



(1 + *)' 



(13) 



Note that in writing this constraint we already combined the one 
for case (III) with the one for case (III.5). It is done by using C\ in 



the lower limit and C2 in the upper limit, such that the combined 
constraint is conservative. From now on we expand the case (III) 
definition to be v ov t < min(v a , v m ) so that it includes case (III. 5). 

Two last pieces of constraining information for case (III) 
are from v a , i.e., v a (f Up ,/3i,R) < vx and v a {f Vv ,fl\,v m ,R) 
> Vopt, where we use the expressions for v a given in Eq. (6); the 
first v a expression is for the situation of v m < v a and the second 
for the situation of v a < v m . The first constraint gives 



Ri4,B c 



x 10 



-1.75-/3 



1 f 2 



y p,19 



( DL,2S \ 

\l + zJ 



(14) 



The v a expression in the second constraint, v a > u opt , contains 
Vm. which can be expressed in terms of f v t and R from the R- 
expression for this spectral case. After substituting for v m , the sec- 
ond constraint gives 



. 1 



x 2.0 



2.4 x 10 P1 



2/3-20! 1 



f 2 

J "apt 



D L 



0- + Z) 



(15) 



Here the constraint v a > v ov t gives a lower limit of R, contrary to 
what is inferred from the conventional relation between v a and R. 
It is because in this subtle occasion v a depends not only on R but 
also on v m , and v m is expressed in terms of the optical flux density 
and R, hence the combined independence of v a is positive. 

The final constraint for case (III) should be the overlapping 
region among those three constraints obtained. 

For case (IV), v m can be obtained from the flux density ratio 
relation by 
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E CM 




Int. 1 (252-282 s) 
X 3 

Int. 2 (402-432 s) 



GRB 050820a 



Int. 3 (515-545 s) 
X 1/3 




10 10'° 10" 10'° 10 
v (Hz) 



1Q 15 10 16 1Q 17 1Q 18 1Q 19 1Q 20 
V (Hz) 



GRB 051 1 1 1 




10 ,J 10'° 10" 10'° 10'" 
V (Hz) 

Figure 2. Observed optical to 7/X-ray broadband spectra for the 4 GRBs in our prompt optical detection sample. The dotted line is a line connecting the 
optical and the lowest energy 7/X-ray data points. For GRB 051111, only one 7-ray data point is shown, superposed with the single-power-law fitted spectral 
index and its confidence range. Adapted from Vestrand et al. (2005), Vestrand et al. (2006), Yost et al. (2007a) and Page et al. (2007), respectively. 



^771,17 



6 x lO 2 " 1 



f, 



p* \ Pi 

U p,W 



01-1/3 



(16) 



whose value will be used to justify the case definition v pt < 
v m < vx- Another constraint from the case definition is 
v a[fv p ,Pi,Vm,R) < v op t. Substituting with the v m expression, 
this gives 



2.4 x 10 



Pi 



fu 



2 £1 
5i 1 , , 2 

P,19 







2/3-2/3! 1 



Dl,2» 



Notice that two contrary constraints, v a > v pt in case (III) and 
v a < v opt in case (IV), give exactly the same constraints on R. 
This is because the v a expression in both cases contains R and v m , 
but in case (III) v m is a strong function of R, oc i? -4 , while in 
case (IV) v m is a function of the flux density ratio only. Thus in 
case (III) the independence is reversed between two sides of the 
inequality relation v a > Vopt ■ 

To summarize, the overall constraints on R are: Eq. (9) for 
case (I), Eq. (1 1) for case (II), Eq. (13-15) for case (III) and Eq. (17) 
for case (IV). In addition, when they are available, the calculated v a 
or v m must satisfy the case definitions. 

If the optical flux density has only an upper limit, the above R- 
constraints must be taken with a conservative point of view wher- 
ever f Vopt is involved. Let f Vopt represent the measured upper 
limit. For case (I), Eq. (9) will give an upper limit for 7?. For case 
(II), Eq. (11) remains. For case (III), Eq. (13) is left with only the 
upper limit of R, Eq. (14) remains and Eq. (15) is useless. For case 
(IV), Eq. (16) gives a lower limit of v m which can be used to justify 
the case definition; Eq. (17) is useless. But we recall that Eq. (17) is 
obtained by substituting the v m expression into the definition con- 
straint v a (f Vp ,f3i,v m ,R) < v op t- Here, instead of using the v m 
expression, we plug in the upper boundary of v m : v m < vx , then 
a new lower limit of 7? is obtained for case (IV): 



R14B, 



i(i±i) 3 >(f)%0, 4 x(0^ 



(18) 



5 GRB DATA SAMPLE 

Now we turn to the real GRB data to which our method devel- 
oped above can be applied. First we construct a small sample of 
GRBs whose prompt optical emission is not only detected but is 
also variable and temporally correlated with the 7-rays. Excluding 
GRB 0803 19B that requires a new spectral component for optical 
emission, we identify four of GRBs in the sample, all belonging to 
the long-duration class. Three of them show complex fluctuations 
in their prompt 7-ray and optical LCs while the fourth is a single, 
smoothly peaked event, so we utilize multiple time intervals for 
each of the three. The emission properties of each time interval are 
listed in Tab. 1. 

GRB 041219A is a very long (T dur ~ 500 s) burst and has 
multiple peaks in 7-ray LC. It has three optical detection intervals; 
the first two are correlated with the first 7-ray peak and the third 
with the second 7-ray peak (Vestrand et al. 2005). 

GRB 050820A is a similar one except that it has denser opti- 
cal temporal coverage. Its optical LC is decomposed into two com- 
ponents : a smooth component with fast rise and power-law de- 
cay, and a strongly variable component superposed on it (Vestrand 
et al. 2006). The smooth component is well accounted for by the 
early afterglow due to the GRB outflow interacting with the ambi- 
ent medium. The variable component is found to correlate with the 
7-ray peaks, suggesting it has the same origin as the 7-rays. In Tab. 
1 the optical emission properties for this burst are for the residual 
optical component after subtracting the smooth component. 

GRB 05 1 1 1 1 has a single FRED (fast rise and exponential de- 
cay) peak in 7-rays lasting ~ 90 s. The first optical observation 
starts at 30 s after the burst trigger when the 7-ray LC began to 
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decay (Yost et al. 2007a). The prompt optical LC (before the 7- 
rays die off) decays more steeply than that of the later optical after- 
glow emission. The prompt optical emission has an excess above 
the back extrapolation of the later optical afterglow component. It 
has a decay slope statistically compatible with that of the 7-ray LC, 
and its flux density is also compatible with the spectral extrapola- 
tion of the 7-rays (Yost et al. 2007a). This is good evidence that the 
prompt optical excess has the same origin as the 7-rays. We use the 
flux density of the excess - not the total - optical emission in the 
first optical observation time interval. 

GRB 061 121 has two separate 7-ray peaks. The last peak was 
caught by XRT, UVOT and ROTSE, and it appears in LCs in all 
bands (Page et al. 2007). We use the emission properties of two 
time intervals of the last peak, one during the rising phase, and the 
other just at the peak. 

We show the SEDs for all time intervals of the 4 GRBs in 
Fig. fusing the data adapted from their original publications. Two 
SEDs show almost no break between the optical and the 7-rays 
(corresponding to the theoretical spectrum case II) or a break very 
close to the optical. In all other cases, at least one break is needed 
between the optical and the 7-rays bands. The break(s) could be v a , 
v m , or both. 

Besides this first data sample for optical detections, we also 
define a second data sample which is composed of those optically 
"dark" GRBs during the prompt phase. This sample is adopted from 
Yost et al. (2007a), who reported the bursts whose prompt phase 
was observed by ROTSE but only upper limits on the optical flux 
were retrieved. Each burst has either a single or multiple time inter- 
vals of ROTSE exposure during the prompt phase. For bursts with 
multiple optical time intervals, we use the interval which has the 
smallest measurement error in the 7-ray flux density f v - usually 
the interval that has the brightest 7-ray flux. The only exception is 
GRB 061222A, for which three time intervals are used. This is be- 
cause all the three intervals are located at the brightest part of the 
7-ray LC. They all have small errors in /„ , and the instantaneous 
7-ray spectral index /?i are available for all three intervals. The 
motivation of selecting this sample is the following: Even if there 
is no direct detection of prompt optical emission, one can speculate 
the existence of a prompt optical emission component that tracking 
the 7-ray LCs, which is the spectral extension of the prompt 7-ray 
spectrum into the optical band. The flux level of this component 
must be fainter than the upper limit set by the ROTSE observations. 
We want to check whether inferred R constraints of this sample is 
consistent with the sample with optical detection, and whether the 
non-detection of optical emission of this sample is due to stronger 
synchrotron self absorption associated with a smaller R. 



6 RESULTS 

We apply the constraints on R derived in Sec. [4] to the first sam- 
ple with prompt optical detections. The results for all four broad- 
band spectral cases are listed in Tab. 1. For most bursts in the 
sample, case (II) can be immediately ruled out because usually 
Popt-x > Pi - Case (IV) is also ruled out for some bursts because 
the calculated v m is > 0.3 keV. For GRB 041219A, case (I) can 
be ruled out for its first time interval because the case definition is 
not satisfied by the calculated v a . Actually this interval is consis- 
tent with case (II), i.e., the optical intensity is consistent with the 
simple power-law extrapolation from the 7-ray spectrum. 

We plot the permitted _R-ranges for each observation time in- 
terval in the sample for all possible spectral cases as floating bars 



in Fig. [3] The observed optical to 7-ray SED restricts v a from be- 
ing much larger than v ovt . Accordingly, the results in Fig. [3] give 
a constraint on the emission site for most time intervals of this 
sample: R ^ a few x 10 14 Ff/^B}/ 4 cm. For two time intervals 
(041219A Int. 3 and 050820A Int. 3) in the sample, some spectral 
cases can be ruled out, thus the i?-constraint can be pinned down 
toR^ (10 14 -10 15 ) r^ 4 B 5 1/4 cm. 

Similar results for the sample with only prompt optical upper 
limits are plotted in Fig. 4. In about half (6/13) of the sample a 
heavy self absorption, i.e. large v a , is needed to account for the op- 
tical deficit, corresponding to the spectral case (I) and (III), which 
implies a constraint of R < 10 15 T^gS^ 4 cm. For the remaining 
half (7/13) of the sample, a spectral break at v m which is below 
ux but is much larger than v ov t alone can give rise to the required 
deficit in optical, while v a can keep being smaller than v op t, cor- 
responding to case (IV). Thus for this half of the sample, we can 
provide no constraint on v a and hence on R. 

Comparing Fig. 4 with Fig. [3] we find that there are always 
overlapping regions between the permitted i?-ranges for the two 
samples. Therefore we can not draw any statistically significant 
distinction between these two samples as regards the constraints 
on their emission sites. 

6.1 The dependence of results on F and B 

Strictly speaking, the constraint on R is dependent on the source LF 
r and the magnetic field strength B. Independent determinations 
of r and B for each GRB in our sample are not easy. So in this 
work we adopt the theoretically anticipated values in the standard 
internal-shock model: Y ~ 300, B « 10 G. In the following we 
will justify these adopted values based on the available information 
of the GRBs in our sample. 

6.1.1 Constraints on T 

Recently Molinari et al. (2007) inferred T w 400 for two GRBs by 
directly observing the deceleration time of the GRB outflow. This 
value is in agreement with what we adopt. In addition, here we 
present some attempts to estimate T for each GRB in our first sam- 
ple (with optical detections) using three independent arguments, 
which suggests that the choice of T ~ 300 is reasonable. 

(1) The variability time scale argument. For the standard sce- 
nario in which the variability time scale is that of the central en- 
gine and that the ejecta form a conical jet with opening angle much 
larger than 1/r, the observed variability time scale should be at 
least the angular spreading time. This gives the constraint 

T>41[R 14 (l + z)/ST] 1/2 , (19) 

where ST is the observed variability time scale, defined as the rising 
or decaying time scale of the pulses in the LC. From the data we 
find that ST is ~ 1 - 15 s for our sample. The T-constraint from this 
argument is rather weak compared with the other two constraints 
described below. 

(2) The photon annihilation opacity argument. A GRB re- 
leases a huge amount of energy (~ 10 53 erg isotropically) in terms 
of high energy photons from a small volume, which are subject to 
the photon-photon annihilation (e.g. Baring & Harding 1997; Lith- 
wick & Sari 2001). Suppose e max is the maximum photon energy 
detected in the burst, the fact that a single power-law or a piece- 
wise power-law spectrum is detected for most GRBs up to e max 
implies that the optical depth of photon-photon pair production r 77 



Table 1. Constraints on the emission site radius for GRBs in which a prompt optical component temporally correlated with the 7-ray emission is 
present. /„ and f3i are the BAT flux density and the BAT spectral index, respectively; they are determined directly from the same time interval as that 
of the optical observation. A Ho =71, = 0.73, £Im= 0.27 universe is assumed. For GRBs without a known redshift, z = 2 is assumed. The "\" 
symbol means the SED is apparently inconsistent with, or the derived v a or u m violates, the spectral case definition, thus the corresponding spectral 
case is ruled out. Spectral case (II) is ruled out for the most of the sample because f3\ is much steeper than (3 op t-x/- r 





time int. 

(s) 


z 


(mJy) 


Up 

(10 19 Hz) 


(mJy) 


Pi 


Reference* 


(I) 


R(10 14 
(II) 


Ti3/4 D l/4 , 
r 300 B 5 Cm ) 

(HI) 


(IV) 


v a (eV) 
(I) 


u m (eV) 
(IV) 




203 - 275 




3 


1 


0.65 


-0.16 




Tfcsl 


> 9.0 


(8.1, 8.7) 


> 8.1 


X 


2.2 


041219A 


288 - 318 


(2) 


10 


1 


3 


-0.32 


1 


18.4 


\ 
\ 


(8.0, 18.4) 


> 8.0 


4 


40 




330 - 403 




4 


1 


2 


-0.64 




14.4 


(3.5, 14.4) 


>-30L 


12 






252 - 282 




7.3 


1.6 


0.56 


-0.3 




15.5 


> 20 


(12, 15.5) 


> 12 


2.4 


4.7 


050820A 


402 - 432 


2.6 


1.3 


1.6 


0.4 


-0.3 


2 


6.5 


\ 
\ 


(2.8, 6.5) 


> 2.8 


4 


42 




515-545 




0.25 


2.6 


0.13 


-0.7 




3.7 


(0.9, 3.7) 




18 


T7G0- 


051111 


32-32 


1.55 


8.2 


1.7 


0.08 


-0.48 


3,4 


18 


> 22 


(14.5, 18) 


> 14.5 


2.3 


3.2 


061121 


56-56 


1.3 


1.8 


1.2 


0.14 


-0.42 


5 


7.8 


\ 
\ 


(4.1,7.8) 


> 4.1 


3.5 


19 


76-76 


7 


1.2 


7 


-0.21 




13 


(4.7, 13) 


> 4.7 


4.3 


97 



* 1: Vestrand et al. (2005); 2: Vestrand et al. (2006); 3: Yost et al. (2007a); 4: Hill et al. (2005); 5: Page et al. (2007). 
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Figure 3. The constraints on the emission radii for 4 GRBs with prompt optical detections. For 3 bursts of the sample, more than one time intervals per burst 
are used. Based on the emission spectral information of individual bursts or time intervals, spectral case II and IV, respectively, are already ruled out for miner 
parts of the sample. See details in Tab. 1. 




Figure 4. The constraints on the emission radii for GRBs without prompt optical detection. GRB 061222A has three time intervals that have information 
available for our calculation. Based on the emission spectral information, spectral case II is already ruled out for the whole sample, and spectral case (IV) is 
ruled out for about half of the sample. For 8 bursts without known redshift, z = 2 is assumed. 



for the photons of energy e max is smaller than unity. r 77 is re- 
lated to the total number of photons and the size of the emission 
region, the latter of which can be expressed in terms of V and 8T 
within the internal shock model. Thus this requirement can impose 
a lower limit on V within the internal shock model. More gener- 
ally the opacity argument can give a constraint in the R — V space 
(Gupta & Zhang 2008; Murase & Ioka 2008; Zhang & Pe'er 2009). 
Our following treatment applies to the internal shock model that is 
commonly discussed in the literature. 

The minimum photon energy at which the photons are able to 
annihilate with photons of energy e ma x is: 



&max \l- ~r Z) 

The power-law form spectrum just below e, nax is described as in 

N v= N ^*AvY a > (2l) 

where e is the detected photon energy, iV(e) - in units of 
[s _1 cm _2 MeV _1 ] - is the number of photons detected per unit 
time per unit area per unit energy at e, iVi is equal to N(e) at e = 1 
MeV, and a is the photon index. 

We followed the formulae of Lithwick & Sari (2001) in their 
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Table 2. The r constraints from the requirement that the photon annihilation optical depth t 77 < 1. A Hq = 71, S7a = 0.73, £Im- 0.27 universe is assumed. 
For the GRB without known redshift, z = 2 is assumed. 



GRB 


time int. 

(s) 


z 


8T 
(s) 


Ni 

(s- 1 cm- 2 MeV- 1 ) 


a 


(MeV) 


&max ,an 

(T| MeV) 


Reference* 


r 




203 - 275 


(2) 


15 


0.15 


1.80 


0.2 


1450 




> 56 


041219A 


288-318 


(2) 


5 


0.56 


1.75 


0.2 


1450 


1 


> 93 




330 - 403 


(2) 


2 


0.10 


2.15 


0.2 


1450 




> 49 




252 - 282 


2.6 


15 


0.30 


2.52 


1.2 


168 




> 53 


050820A 


402 - 432 


2.6 


15 


0.063 


2.00 


0.5 


403 


2 


> 53 




515-545 


2.6 


10 


0.015 


1.96 


0.25 


806 




> 41 


051111 


32-32 


1.55 


8 


0.035 


1.48 


0.15 


2680 


3 


> 65 


061121 


76-76 


1.3 


1 


2.7 


< 2.9 


1.4 


359 


4 


> 53 



* 1: Vestrand et al. (2005); 2: Cenko et al. (2006); 3: Yost et al. (2007a); 4: Page et al. (2007). 



Table 3. The T constraints derived from the deceleration time tdec constraints for ISM and Wind medium, respectively. A Ho =71, f2^ = 0.73, Om= 0-27 
universe is assumed. For the GRB without known redshift, z = 2 is assumed. For GRB 041219A the inferred idec-constraint is very loose, partly because of 
a lack of X-ray afterglow observation, and also because the early infrared light curve (( ^ 6 X 10 3 s) is highly variable, possibly of internal shock origin, 
which makes it difficult to infer tdec to be these earlier times. In the cited reference for GRB 061 121, only i? 7i ,: so is given, without giving the 7-ray fiuence. 



™„ tdec 7-ray fiuence E lAao T 

GRB z Reference , /a , ,„ , /A 

(s) (IO- 5 erg cm- 2 ) (10 53 erg) (*7 7 ,o.2 ™ /8 ) fa 7 ,o.S >. A * > 

041219A (2) < 3 x 10 4 15.5 15.2 1 > 73 > 33 

050820A 2.6 < 500 5.3 8.3 2, 3 > 337 > 81 

051111 1.55 < 100 0.39 0.24 4 > 349 > 46 

061121 1.3 < 200 - 2.8 5 > 352 > 70 



* 1: Vestrand et al. (2005); 2: Vestrand et al. (2006); 3: Cenko et al. (2006); 4: Yost et al. (2007a); 5: Page et al. (2007). 



Limit A case to calculate the lower limit of T due to r 77 < 1 for 
our optical detection sample. The results, as well as the observa- 
tional properties that are used, are summarized in Tab. [2] Note that 
for all bursts in the sample, e max ,an 3> e ma x for reasonable values 
of T (e.g., ~ 10 2 ).Thus our calculation has implicitly assumed that 
the power-law form spectrum detected below e max would actually 
extend well beyond e max and up to e m ax,an, which is usually be- 
low the detector's limited bandpass. More rigorous calculations re- 
quire the knowledge of the spectral shape below peak energy of 
the spectrum (Gupta & Zhang 2008). However, the optical depth is 
much more sensitive to T than the spectral indices, and the derived 
F constraint is not significantly modified by performing the more 
rigorous treatment. The above treatment is adequate to serve our 
purpose. 

(3) The deceleration time constraint. In the well established 
external shock model for GRB afterglows (e.g., Meszaros & Rees 
1997; Sari, Piran & Narayan 1998; Chevalier & Li 2000; see Pi- 
ran 2005 for a review), the onset of afterglow marks the time, as 
known as the deceleration time tdec, when one half of the total ki- 
netic energy of the GRB outflow Ek is deposited to the shocked 
circumburst medium. The afterglow light curve should rise before 
tdec and decay in a power law with time after tdec (e.g., Moli- 
nari et al. 2007). The deceleration time tdec is determined by Ek, 
outflow bulk LF T and the density of the circumburst medium, 
either a constant density medium (ISM) or a wind-like medium 
lp(R) = A x R- 2 ]: 



tde 



( 3Ek -Z X /3 (l I z) 

\64ivnm p c t 'r 11 J \ L ^ Z ) 
£fc.,.,„(l + z) 



for ISM, 
for Wind, 



(22) 



where n is the proton number density of the ISM medium, and A = 
5 x 10 11 A* g cm -1 is the Wind medium density normalization 
parameter. The isotropic equivalent kinetic energy Ek,iso can be 
related to the isotropic energy release in 7-ray radiation E^.iso by 
an energy conversion efficiency factor ?7 7 = E lt i 30 / Ek,iso- 

The afterglow observations (X-rays and optical) for our sam- 
ple show either a power law decay starting from the earliest obser- 
vation interval or, a long-lasting shallow decay followed by a nor- 
mal power law decay. It suggests that the afterglow onsets should 
be earlier than the start of the single power law decay or the start 
of the shallow decay. One can use the first observation data in the 
decaying afterglow phase to constraint the deceleration time to be 
earlier than the observational epoch (e.g., Zhang et al. 2006). From 
the data, we find the deceleration time tdec < 3 x 10 4 s for GRB 
041219a, and t dec < 100 - 500 s for the other 3 GRBs (050820a, 
051111, 061121). Assuming Tft = 0.2, n = 1 cm" 3 and A, = 1, 
we find: for ISM, V > 73 for 041219a, V £ 350 for the other 3 
GRBs; for Wind, T > 33 for 041219a, F > 50 for the other 3. 
The results are summarized in Tab. [3] We have checked of the com- 
pliance with closure relationships predicted by the external shock 
models for these bursts during the afterglow phase. It turns out that 
two (GRB 050820A and 061121) out of the 4 GRBs are consis- 
tent with and in favour of the ISM environment scenario, while the 
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other two are consistent with both scenarios and can not discrimi- 
nate between them. 

To summarize, with the available data, one can only constrain 
but cannot determine T of the GRBs in our sample. On the other 
hand, all three constraints derived from the data are consistent with 
r = 300 adopted in our calculations. In particular, the T-constraint 
derived from the argument (3), which is the most stringent one 
among the three arguments, indicate that the assumed value of 
r = 300 is reasonable. 



6.1.2 Constraints on B 

The B value in the emission region is a function of R. There are 
two possible origins of the magnetic field in the emission region. 
The first component is the global magnetic field entrained by the 
ejecta from the central engine. Let's assume B ~ 10 14 G at the 
central engine, a typical value for a fast rotating magnetar or a fast 
rotating black hole accretion disk system - the two most plausible 
GRB central engine candidates. The B value drops as R~ 2 with 
R up to the light cylinder, and then drops as R^ 1 thereafter (Gol- 
dreich & Julian 1969). Given that the central rotating source has a 
radius of R t ~ 10 6 cm and a rotation period of P ~ 1 ms, at a 
radius R ~ 10 14 cm, the field that is carried within the outflow has 
a strength of B « 2 x 10 s B^uRl^P^l G. 

The second B component is a random field generated in-situ 
in the emission region, likely in a relativistic shock via the Weibel 
instability (Medvedev & Loeb 1999). This random field also fol- 
lows the same independence and is of the same order as the engine- 
related B component if e_g - the ratio of the post-shock magnetic 
energy density to the total energy density - is not too small (Zhang 

6 Meszaros 2002). 

One can also briefly estimate the local field strength by relat- 
ing the GRB 7-ray peak photon energy, typically ~ 0.1 MeV, with 
the synchrotron characteristic frequency, v v = - njr^ ■ where 

7 is the typical energy of electrons. In the internal- shock model, 
7 = e e f(p)(m p /m e )9 p , where e e is the ratio of the electron en- 
ergy density over the total thermal energy density in the post-shock 
fluid, f(p) — (p — 2)/(p — 1) and p is the electron energy spectral 
index. The parameter 9 P is the fractional energy gain of a proton 
passing the shock which depends only on the relative LF between 
the fluids downstream and upstream. For internal shocks, 8 P is not 
dependent on the shell bulk LF V and is of the order of unity. For 
e e = 0.3 and p — 3, we have 7 ~ 300. Thus the required B value 
can be estimated as B w 5 x W 5 ^Joo"f30oi^') G. 

The three crude estimates are marginally consistent with each 
other. We have taken B = 10 5 £>5 G as the typical value throughout 
the text. Of course a large uncertainty exists due to our lack of 
understanding on the field properties, but it is reconciled by the 
very weak dependence of R on B (1/4 power). 



6.2 Comparison with results from an alternative modelling 
approach 

Kumar & McMahon (2008) developed a general method of mod- 
elling GRB's 7-ray emission properties. Their method considers 
the synchrotron and the SSC emission, respectively, as the radia- 
tion mechanism, takes into account the radiative cooling of elec- 
trons, and uses observed emission properties (such as the peak flux 
density and the pulse duration) to search for the allowed space of 
the model parameters such as R and T. Here we also apply their 



method for the synchrotron case to our optical detection sample us- 
ing their code, and compare the results with ours. We add a new 
constraint into the module that controls the allowed model parame- 
ter space, which is that the optical flux density calculated from the 
model has to match the observed one within a ±50% range. 

To use this detailed modelling method, we have to specify 
which standard synchrotron spectral regime a GRB emission inter- 
val is in. The spectral indices (J3i) of the optical detection sample in 
Tab. 1 have a variety of values around -1/2, based on which we clas- 
sify the sample into 4 categories and apply the detailed modelling 
method accordingly. ( 1 ) For those time intervals that are most prob- 
ably consistent with the f3i = -1/2 regime (051111 and 061121 
int. 1), this method gives R ss 10 14 — 10 16 cm. (2) For those possi- 
bly consistent with both (3\ = — 1/2 and (3i = — (p — 1 ) /2 where 
p > 2 (041219A int. 3 and 050820A int. 3), the method in the 
/?i = -1/2 regime gives R w 10 14 - 10 15 cm, while in the fix = 
— (p — l)/2 regime it gives no allowed i?-space — but if we relax 
the j v t constraint, it gives R « 10 17 — 10 18 cm. (3) For those 
possibly consistent with both (3i = —1/2 and (3\ = — (p — l)/2 
where 1 < p < 2 (041219A int. 2, 050820A int. 1 and 2), the 
method in the f3i = -1/2 regime gives R w 10 14 - 10 15 cm; 
however, the code provided by the authors is not applicable when 
p < 2. (4) The last category are those inconsistent with j3i = — 1/2 
but probably consistent with f3i = — (p — l)/2 if 1 < p < 2 
(041219A int. 1 and 061121 int. 2) for which the code is not appli- 
cable. 

Overall, we find that for the major part of the optical de- 
tection sample where the detailed modelling method (Kumar & 
McMahon 2008) is applicable the allowed spaces for R from this 
method are about 10 14 - 10 16 cm. This is approximately con- 
sistent with the findings from our approach that R is ^ (a few 
x 10 14 - 10 15 ) T%^Bl /A cm for most of the intervals in the sam- 
ple, and 10 14 r^Bs 74 cm < R < 10 15 T%qBI /a cm for the 
remaining two intervals in the sample. 



7 CONCLUSION AND DISCUSSIONS 

Based on the assumption that the prompt optical and 7-ray emis- 
sions belong to the same synchrotron continuum of a group of 
hot electrons, we make constraints on the location of the prompt 
emission site for a sample of GRBs whose prompt optical emis- 
sion is detected to temporally tracking the 7-ray LCs, by de- 
termining the location of v a in their SED. Our analysis shows 
that for most of the intervals in this sample the distance of the 
prompt emission site from the explosion centre R is ^ (a few 
x 10 14 — 10 15 )r3Q 4 _Bg^ 4 cm, and for the remaining two intervals, 
the emission site is (10 14 — lO^r^i?,^ 4 cm away from the ex- 
plosion centre. 

The dependence of the distance constraint on the GRB outflow 
LF r is not negligible. On the other hand, various indirect observa- 
tional constraints on V point to T ^ 300 (e.g. Molinari et al. 2007; 
and even T > 600 for GRB 080916C, Abdo et al. 2009). The de- 
rived observational constraints on T for bursts of this sample (Sec. 
16. Il l are consistent with such an inference. In our paper, we take 
F = 300 as a typical value. A higher F would only make the above 
distance constraint even larger. Our knowledge of the local mag- 
netic field strength B is less certain, although the 7?-dependence of 
B is weak. Several crude estimates of B based on the synchrotron 
radiation mechanism for GRB prompt emission suggest B ~ 10 5 
G at R ~ 10 14 cm. This typical value has been adopted in our 
calculations. 
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The ^-constraint we obtained is inconsistent with the photo- 
spheric emission model in which the prompt emission arises at the 
photosphere radius of 10 11 — 10 12 cm (Rees & Meszaros 2005; 
Ryde et al. 2006; Thompson et al. 2007). This result alone can not 
discriminate between the fireball internal shock model and the mag- 
netic outflow model. By comprehensive modeling the GRB prompt 
7-rays and early X-rays, Kumar et al. (2007) concluded a prompt 
emission site of 7? ~ f 15 — 10 16 cm, which is supported by their 
further general modeling of the 7-ray emission properties (Ku- 
mar & McMahon 2008). A large R is derived for GRB 080916C 
through thepair opacity constraint (Abdo et al. 2009; Zhang & 
Pe'er 2009D and for GRB 080319B through the synchrotron self- 
absorption constraint (Racusin et al. 2008) and the SSC scattering 
optical depth constraint (Kumar & Panaitescu 2008). In summary, 
a large prompt emission distance from the central engine seems to 
be supported by three independent approaches, respectively, i.e. in 
7/X-rays (e.g., Lazzati & Begelman 2005; Lyutikov 2006; Kumar 
et al. 2007), in GeV 7-rays (Abdo et al. 2009; Zhang & Pe'er 2009), 
and in the optical band (this paper). 

We have also studied a sample of GRBs with prompt opti- 
cal non-detections. Applying the same technique, we do not find 
any inconsistency between their R constraints and those of the op- 
tical detection sample. This result is inherited from the findings 
by Yost et al. (2007b) that no distinction in distributions of /3i 
and /3 op i_ 7 can be drawn between the optically dark GRBs and 
the GRBs with optical detections. However, this is only because 
the currently limited instrumental sensitivity prevents a distinction 
from being drawn. Deeper observations in the future in the optical 
band would provide further information regarding whether the op- 
tical deficit is due to a heavier synchrotron self-absorption in these 
bursts. 

7.1 Multi-color information for the low-energy spectrum 

If multi-color photometry near the optical band exists for the same 
time interval during the prompt phase, it would provide the local 
spectral index near the optical band (provided that the extinction 
correction is properly made). This would be helpful to identify the 
spectral case the data satisfy. For example, the spectral indices near 
optical differ by A/3 = 5/3 between the cases (III) and (IV), and by 
A/3 =1/2 between the cases (I) and (III). Unfortunately, this kind 
of observational information is unavailable for all the time intervals 
of the optical detection sample we have considered in Tab. 1 . The 
hope is that future multi-band prompt optical detections may be 
able to break the spectral case degeneracy and to tighten the 7?- 
constraint. 

7.2 Limitations of the method 

Our method is based on the assumption that the optical emission 
is emitted from the same group of electrons that produce the 7- 
rays in the same site via the same synchrotron radiation. This one 
component assumption has some supports (see Sec. 3) but certainly 
is not conclusive. 

There are three other scenarios that have been discussed in 
the literature (mostly motivated by interpreting GRB 0803 19B). 
The first one invokes two emission zones for optical and 7-rays. 

2 Although both work obtained large R, the inference of R in Abdo et al. 
(2009) is based specifically on the internal shock model while Zhang & 
Pe'er (2009) gave a more model-independent constraint on R. 



For example, in the internal-shock-model based residual collision 
scenario proposed by Li & Waxman (2008), the shells with high 
LF contrast in a GRB outflow collide first and merge at smaller 
radii, producing the 7/X-rays. Later those merged shells with low 
LF contrasts would collide mildly at later times and larger radii, 
giving rise to optical emission. Alternatively, if the outflow is neu- 
tron rich, the proton shells tend to collide at smaller radii to power 
7/X-rays, while the free neutrons only decay at large radii and the 
decay products would be collided by later injected faster proton 
shells and power optical emission at larger radii (Fan et al. 2009). 
In both scenarios, it is expected that the observed optical pulse peak 
emission time is delayed by R opt /2T 2 c with respect to the 7-ray 
pulse peak emission time, which may be in principle tested if the 
data quality is high. These models however do not naturally predict 
a smooth extension of 7-ray spectrum to the optical band without 
distinct spectral features. Although they cannot be ruled out by the 
data in our sample, they are more complicated than our one-zone 
model. 

Secondly, our analysis is not applicable in the following sce- 
nario (Yu et al. 2009): a pair of shocks (reverse shock and forward 
shock) arise when two shells collide in the internal-shock model; 
different populations of electrons are accelerated in each of the two 
shocks and the two populations have different typical electron en- 
ergies and different shock-generated magnetic strengths. The char- 
acteristic synchrotron frequencies are different - the forward shock 
produces the optical emission, while the reverse shock produces the 
7-ray emission. These two emissions are two spectrally indepen- 
dent components but are temporally correlated because the heating 
of the two electron populations arises from the same dynamical pro- 
cess. Although this model may interpret the peculiar SED shape of 
GRB 0803 19B whose prompt optical flux density exceeds the ex- 
trapolation from the 7-ray spectrum by 4 orders of magnitude, it is 
unclear whether it can work properly for the bursts in our optical 
detection sample. 

Finally, the synchrotron + SSC scenario (Kumar & Panaitescu 
2008; Racusin et al. 2008) has been proposed to interpret GRB 
0803 19B. For our sample, there is no need to introduce a second 
distinct spectral component. We consider only pure synchrotron ra- 
diation in deriving constraints on R, although we have given the v a 
estimation for SSC radiation whose value is similar to u a for syn- 
chrotron. In principle, it is possible that the observed emission from 
optical to 7-rays is dominated by SSC. If this is the case, our ap- 
proach of constraining R by calculating v a may give results some- 
what different from the synchrotron case for the same sample of 
bursts. We did not carry out the analysis for the SSC case because 
the detailed shape of the SSC spectra is much more complicated 
than the synchrotron one. Note that our approach assumes that the 
optical and the 7-rays are from the same group of electrons due to 
the same radiation mechanism. In the SSC scenario this approach 
is applicable only if the SSC component dominates a large spectral 
band from the 7-ray down to the optical and the synchrotron com- 
ponent has to lie well below the optical band. This is usually not 
expected in the SSC models (e.g. Piran, Sari & Zou 2009). In any 
case, if SSC is involved in interpreting any part of the spectrum in 
our sample, then our analysis is no longer applicable. 
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APPENDIX: DERIVATION OF THE SELF-ABSORPTION FREQUENCY 

In this appendix, we provide a rigorous derivation of the blackbody equivalence equation (Eq. 3) that we use to calculate the self-absorption 
frequency (Eq. 3). The derivation is carried in the GRB ejecta comoving frame in which the relevant quantities are marked with the prime 
sign. 

The self-absorption frequency v' a is defined as T{y' a ) — 1, where r(v' a ) is the optical depth due to the self-absorption at v' — v' a . The 
optical depth t(v') = J a' v ,ds' deceases with the frequency, where a' u , is the self-absorption coefficient [cm -1 Hz -1 ] and the integral is 
over the line-of-sight width of the emitting source. The integral can be calculated directly only when we have the exact information on the 
number density of the emitting particles and its distribution over the length, which is not easy. Instead of directly calculating the integral 
J acids', we turn to derive the emission coefficient j' u , [erg s -1 cm -3 Hz -1 sr -1 ] and then express the integral of ot v , over width into the 
integral of j' v , over width, the latter is just the specific intensity at the source surface which is directly observable. 

The synchrotron radiation spectrum, or specific radiation power, of an electron with LF 7 gyrating in a magnetic field B with a pitch 
angle a is 

P iy ,l) = F ( u v ch), (A-l) 

where e is the electron charge and 

, ZeB sin Q7 2 

«4 = — j (A-2) 

is the characteristic photon frequency of the electron. The function 

F{x) = x / K 5/3 (A-3) 



has the asymptotic form 

f ^ , A x/2) 1 ' 3 ~2.15x 1/3 , iix < 1, 

F x) ~ <^ V3T(l/3) w 1 ' ^ ' 

where T(l/3) is the gamma function of argument 1/3, and it reaches the maximum F max (x) ~ 0.92 at x ~ 0.29. One integral property of 
the function F(x) is 

2 M + 1 



x^F(x)dx - 



i) 



1* + 



where T(y) is the gamma function of argument y. We will use this property later. 

The self-absorption coefficient for any radiation mechanism is given (Rybicki & Lightman 1979) by 



8nm e i/' 2 J ' d'y 



AT( 7 ) 



7 2 



(A-6) 



where P'(v', 7) is the single electron's specific radiation power, and N^d'y is the number density of electrons with energy in the interval 
from 7 to (7 + d'y). In the case of GRB, N(-y) has a two-power-law form and was described in Sec. 2 of the paper. We rewrite it here as 

m = \ C^t-^, if7>>, (A ' 7) 

where C-, is the normalization constant. Note that this distribution set-up is phenomenologically based on the two-power-law shape of the high 
energy radiation spectrum observed in GRBs; in the context of some specific particle acceleration scenario, e.g., shock acceleration, where 
two characteristic electron energies, i.e., the injection energy 7i and the cooling energy y c , are involved, there will be 7 m = min(7i, 7c ) and 
7 P = max(7i, 7p ). Thus d[N( , y)/'y 2 ]/d'y = (— p — 2)N(j) where p could be either pi or P2 depending on the location of 7. For the 
synchrotron radiation (Eg. I A- It the self-absorption coefficient would be 



%/3e 3 B sin aC 7 
87rraic 2 i/ 2 



rip roo 
(pi+2) / F( 2; )7- (pi + 1) d7+(p2 + 2) / F(s) 7 - (P2+I) d 7 

Am J lp 



(A-8) 



where a; = x(j) = v' /v' ch {pf) = (47rm e c^')/(3eB sin Q7 2 ). 

Then we consider two different locations of v': v' < u' m and v' m < v' < v' v , respectively. If v' < v' m , then F(x) falls in the oc x 1//3 
asymptotic regime. One can transform the integral in Eq. dA-8b for the variable 7 into the integral for the variable x. Notice the contribution 
from the second integral part in Eq. dA-8b is unimportant, as long as 7,™ <§C 7 P and 1/3 < pi < p2- Thus it gives 

/ = 1 (Pi + 2) e 3 B sin aC 1 / 47rm e e \ _( Pl+2 /3) ,-5/3 , . q , 

a "' 2 4 / 3 r(l/3) (pi +2/3) m 2 c 2 \3eBsma) 7m ' ( ' 

If v' m < v' < v' p , then x(jy m ) «C 1 and aj( 7p ) ^> 1. After the transformation of the variable 7 into the variable x, the first integral 
part in Eq. JA-8t is in effect integrating over the 2>range from i( 7 m) ~ to i( 7p ) ~ 00, thus we can use Eq. iA-5l to calculate it. For the 
second integral part of Eq. dA-8K f(x) oc x 1 ^ 3 , but its contribution is unimportant as long as 7to <C 7p and 1/3 < p\ < p2- Therefore it 
gives 
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a > v , = /^V^ 2 in (pi+2)/2r (3 P1 + 22N ( 3pi+2^ ^ (P1+4)/2 

87rmlc 2 V27rm e c/ ; V 12 / V 12 / 

To calculate i/„ from T(v' a ) = J a„, ds' = 1, one has to know the exact information about N(j m ) and its instantaneous distribution 
over the width of the emitting source along the line of sight, which are always subject to uncertainties. Nevertheless, in their attempts to 
calculate u' a , some authors have calculated iV(7 m ) by assuming all electrons swept up by the shock are accelerated to relativistic energies 
either in the blastwave model for GRB afterglows (Granot, Piran & Sari 1999; Panaitescu & Kumar 2000; Pe'er & Waxman 2004) or in 
the internal-shock model for prompt emissions (Li & Waxman 2008). We warn that the reality in nature may be that not all but only a 
small fraction of the electrons encountered by the shock can be heated to relativistic energies and radiate, as was suggested by Bykov & 
Meszaros (1996) and Daigne & Mochkovitch (1998) (also see Kumar & McMahon (2008) for an idea of repeated acceleration of a group of 
electrons), and this will introduce the biggest uncertainty to N(-y m ) hence to this "conventional" approach of calculating v' a . In general, this 
conventional approach over-estimates the number of emitting (and absorbing) electrons and, hence, over-estimates v a - In the literature, it is 
usually suggested that v a is slightly below the X-ray band. According to our corrected calculation, v a is typically lower and can extend to 
close to the optical band in a wide parameter range. 

Our new approach here is to express J acids' in terms of J j' v ,ds', both of which contain the term N(j m ) but the latter one is directly 
observable - it is just the specific intensity at the source surface. Therefore the new approach can avoid the uncertainties associated with 
partial acceleration and inhomogeneity over the source radial width. 

Let us calculate the emission coefficient j' v , . By definition, 
rip r°° 

iirjl'= P'(u', 7 )N{~f)dj+ PV, 7 )iV(7)d 7 . (A-ll) 
Am J-y P 

Following the same procedure of calculating ol v , , the integration gives 



1 e 3 BsmaC T (1/3-pi) ,1/3 -r / , 

2 1 /3r(l/3)(p 1 -l/3) ™~7^ \3eBsina) 1™ ' u " ^ "»> 



j' f — J ■2. + fL\L/-A)(pi-l/-i) m B c- \c.ens,Liat/ ■■■■ (A 12) 

The ratio of j' v , over a' u i , also called the source function, is 
& = ^ = ^ (^^^ r(^2)r(^i) (A-13) 

which does not have dependence on N(-y m ). It shows that, for synchrotron radiation, the power-law index of the optical thick (to the self 
absorption) part of the emergent spectrum below v' m is 2, while the power-law index of the optical thick spectrum above v' m is 5/2. 

From the definition of the self-absorption frequency J a u i ds' = 1, we have J (j' v , / S' v i )ds' = 1. Since S' u i does not depend on 
N(fm) and its distribution over the source width, it can be taken out of the integral. Thus we have S' v , = J j' v , ds' = F' v , , where F' v , is the 
specific flux at the source surface in the asymptotic optically thin regime at v' a . Rewriting the expression for S' u , (Eq. IA-131 > at v' = v' a and 
using the photon frequency vs. electron energy relation v' ch (~j) for synchrotron radiation, we get 

max(7 m ,7 a ) x 2m e ^ 2 = F^Cfjn), (A-14) 
where 7 a is the energy of the electron whose characteristic photon frequency is v' a , and the correction factor 

r c x (p x ) = (pi+ ; l ) | p 2 i / 3 1/3) , if < < v'm, 

C<Pi)={ r , , ft, , 1 , r(S£i+H)r(%a) , (A-15) 

y Ga(pi) = V2(pi + l) r ^ 3pl+19 j r ^ 3pl -i^ lf v m <Va< "p. 

If we assume a temperature T' = max( 7m , y a )m e c 2 /k then the last equation is 

2kT'^ = FUCfr), (A-16) 
where for practical uses the correction factor C(pi) is changed to C(J3\) using the relation f3\ = —(pi — l)/2, and so 

C(/3l)= | C 2 (/3x) = 2^(l-/3x)^S|^|S} J iiv' m <v' a<V ' p . ■ (A " l?) 

In the samples presented in the main body of the paper, f3i is among -1.4 to 0, so the ranges for the correction factor are Cx(/3x) = (1-2, 4.5) 
and C2(/32) = (1.2, 7.0). Therefore Eq. J A- 16b shows that, within a factor of a few, at v' a the un-absorbed source surface flux density is 
equal to the flux density of the Rayleigh- Jeans part of the blackbody spectrum with a temperature corresponding to the lowest energy of 
those electrons that are barely affected by the self absorption. This equation is used to calculate v' a in the main body of the paper where the 
correction factor C(j3i) is taken into account. 



